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A search for the non- conservation of lepton flavor in the decay r — > e ± j has been performed 
with 2.07 x 10 s e + e~ —* t + t~ events collected by the BABAR detector at the PEP-II storage ring 
at a center-of-mass energy near 10.58 GeV. We find no evidence for a signal and set an upper limit 
on the branching ratio of B(t — * 67) < 1.1 x 10 -7 at 90% confidence level. 

PACS numbers: 13.35.Dx, 11.30.Hv, 14.60.Fg 



Lepton flavor conservation differs from other conserva- 
tion laws in the Standard Model (SM) because it is not 
associated with an underlying conserved current symme- 
try. Consequently, new theories attempting to describe 
nature beyond the SM often include lepton flavor violat- 
ing processes such as the neutrino-less decay of a /1 or r 
lepton, which have long been identified as unambiguous 
signatures of new physics. If no specific theoretical model 
is assumed, any or all of the /i — > e7, t — > /i 7 and r — > ej 
decays can be expected to be observed, and therefore in- 
dependent searches for each of these modes are required. 
Some theoretical models Q, respecting the current lim- 
its on BO+ -► e+ 7 ) and -> ^±7) Q in fact 
allow r — > e^7 decays to occur up to the existing ex- 
perimental bound |5j. 

A significant improvement on this — ► e^7 limit is 
presented here using data recorded by the BABAR detec- 
tor at the SLAC PEP-II asymmetric-energy e + e _ storage 
ring. The data sample consists of an integrated luminos- 
ity of L= 210.6 fb _1 recorded at a center-of-mass (cm.) 
energy (y/s) of ^/s — 10.58 GeV, and 21.6 fb _1 recorded 
at y/s = 10.54 GeV. With an average cross section of 
tTe+e-^r+r- = (0.89±0.02) nb @ as determined using 
the KK2F Monte Carlo (MC) generator @], this corre- 
sponds to a data sample of 2.07 x 10 8 r-pair events. 

The BABAR detector is described in detail in Ref. ||. 
Charged particles are reconstructed as tracks with a 5- 
layer silicon vertex tracker and a 40-layer drift chamber 
(DCH) inside a 1.5 T solenoidal magnet. An electro- 
magnetic calorimeter (EMC) consisting of 6580 CsI(Tl) 
crystals is used to identify electrons and photons. A ring- 
imaging Cherenkov detector (DIRC) is used to identify 
charged hadrons and provides additional electron identi- 
fication information. 

The signature of the signal process is the presence of 
an isolated ej pair having an invariant mass consistent 
with that of the r (1.777 GeV/ c 2 @|) and a total energy 
(E ei ) equal to y/s/2 in the cm. frame, along with other 
particles in the event with properties consistent with a 
SM r decay. Such events are simulated with higher-order 
radiative corrections using the KK2F MC generator Ml 
where one r decays into e7 according to phase space [lOj , 



while the other r decays according to measured branch- 
ing ratios [llj simulated with the TAU0LA MC genera- 
tor 0, E| . The detector response is simulated with the 
GEANT4 package |l4j. The simulated events for signal as 
well as SM background processes [IEEE EE HQ are 
then reconstructed in the same manner as data. The MC 
backgrounds are used to optimize the selection criteria 
and study systematic errors in the efficiency estimates, 
but not for the estimation of the final background rate, 
which relies solely on data. For the background from 
Bhabha events, we do not rely upon MC predictions be- 
cause the large Bhabha cross section makes generation of 
a sufficiently large MC sample impractical. 

Events with zero total charge and with two or four 
well-reconstructed tracks inconsistent with coming from 
a photon conversion are selected. The event is divided 
into hemispheres by the plane perpendicular to the thrust 
axis. The thrust axis, which characterizes the direction of 
maximum energy flow in the cm. frame of the event |l8j . 
is calculated using all observed charged and neutral par- 
ticles. 

The signal-side hemisphere is required to contain at 
least one 7 with a cm. energy greater than 500 MeV, and 
one track identified as an electron. The electron identifi- 
cation uses DCH, EMC and DIRC information, including 
a requirement that the E/p ratio (the energy deposited in 
the EMC by the charged particle divided by its momen- 
tum as measured in the DCH) lies between 0.89 and 1.2. 
The electron candidate is required to lie within the fidu- 
cial acceptance of the EMC and to have a momentum 
greater than 500MeV/c. These criteria yield a it mis- 
identification rate of less than 0.3%. The efficiency for 
correctly identifying reconstructed tracks in the fiducial 
volume as electrons in r ± — * e ± 7 MC events is greater 
than 91%. For events with more than one signal-side 7 
candidate, we choose the 7 which gives the mass of the ej 
system closest to the r mass. This provides the correct 
pairing for 99.9% of selected signal MC events. 

The resolution of the cy mass is improved by as- 
signing the point of closest approach of the e track to 
the e + e~ collision axis as the origin of the 7 candidate 
and by using a kinematic fit with E ei constrained to 
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y/s/2. The resulting energy-constrained mass (toec) and 
AE = E ei — \fs/2 are independent variables apart from 
small correlations arising from initial and final state ra- 
diation. The mean and standard deviation of the Toec 
and AE distributions for reconstructed MC signal events 
are: (m EC ) = 1777MeV/c 2 , a(m EC ) = 9MeV/c 2 , (AE) 
= -15MeV, a(AE) = 51MeV, where the shift in (AE) 
comes from photon energy reconstruction effects. To 
minimize possible biases, we perform a blind analysis by 
excluding all events in the data within a ±3tr rectangu- 
lar box centered on (toec) and (AE) until all optimiza- 
tion and systematic studies of the selection criteria have 
been completed. We optimize the selection to obtain the 
smallest expected upper limit in a background-only hy- 
pothesis for observing events inside a ±2<r rectangular 
box signal box defined by: \AE - (AE)\ < 2a(AE) and 
|toec — ™t| < 2<7 (mEc ) j as shown in Figure 

The dominant backgrounds arise from Bhabha and 
e + e~ — > t + t~ (with a r — > evv decay) processes with an 
energetic 7 from initial or final state radiation or from 
t — * evvj decays. Backgrounds arising from radiation 
are reduced by requiring that the total cm. energy of all 
non-signal 7 candidates in the signal-side hemisphere be 
less than 200 MeV. To suppress non-r backgrounds with 
significant radiation along the beam directions, the polar 
angle (9 miss ) of the missing momentum associated with 
the neutrino(s) in the event is required to lie within the 
detector acceptance (—0.76 < cos6 m i SS < 0.92). 

The tag-side hemisphere, defined to be that opposite 
to the signal-side hemisphere, is expected to contain a 
SM r decay characterized by the presence of one or three 
charged particles and missing momentum due to unob- 
served neutrino(s). Taking the direction of the tag-side r 
to be opposite the signal ej candidate, we use all tracks 
and 7 candidates in the tag-side hemisphere to calculate 
the invariant mass squared of the tag-side missing mo- 
mentum (to 2 ,), which peaks around zero for the signal. 
To reduce backgrounds from radiative e + e~ — ► t + t~ 
processes, we require to 2 > —0.25 GeV 2 /c 4 . 

The component of the missing momentum of the event 
transverse to the collision axis scaled to the beam en- 
ergy (2 x Pmiss/V^) 1S expected to be large for signal 
and e + e~ — ► t + t~ events, but small for Bhabha and 
2-photon events. We exploit an observed correlation be- 
tween to 2 and (2 x p^iss/V^) m the non-r backgrounds 
to significantly suppress them. We require the following: 
(TO 2 /1.8GeV 2 /c 4 ) - ln(2 x p T mls J J7s) / 2.0 < 1, the high- 
est cm. momentum track on the tag-side hemisphere to 
be inconsistent with being an electron, including require- 
ments that E/p be less than 0.5 and that the momentum 
be greater than 500 MeV/c, and the tag-side hemisphere 
to have a total cm. momentum of all charged and neutral 
particles less than 4.75 GeV/c. 

Backgrounds from e + e _ — > qq processes are further 
reduced by requiring the total invariant mass of particles 
in the tag-side hemisphere to be less than 1.8GeV/c 2 . 
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FIG. 1: r?iEC vs. A_E distribution of data (dots) and shaded 
region containing 50% of the selected signal MC events inside 
the Grand Signal Box, as defined in the text. The boundary 
of the ±2<t signal box is also shown. 



After this selection, 8.9% of the total generated MC 
signal events survive within a Grand Signal Box (GSB) 
region defined as follows: toec € [1.5, 2.0] GeV/c 2 , AE S 
[—1.0,0.5] GeV. The data distribution of toec and AE 
inside the GSB is plotted as dots in Figure ^ along with 
a shaded region containing 50% of the selected signal 
MC events shown for illustrative purposes. The GSB 
excluding the ±3ct blind region contains 1110 data events, 
while the luminosity-normalized sum of the non-Bhabha 
MC backgrounds yield 1045 events. Of these MC events, 
99.8% are e + e~ — > t + t~ events, 99.9% of which have 
t — > evv decays on the signal-side. 

The (5.9±3.7)% difference between the number of data 
and r-pair dominated MC events indicates that the 
Bhabha background level in the GSB is low. However, in 
the more restrictive \ AE - (AE) \ < 2a(AE) region, the 
Bhabha background is expected to contribute a substan- 
tially higher background fraction because of the greater 
likelihood of a Bhabha than a r-pair event to have a 
hemisphere containing the full beam energy. This resid- 
ual Bhabha contamination is studied using data distri- 
butions of the deviation (AE 7 ) of the measured photon 
cm. energy from the corresponding prediction assuming 



a fully contained e + e 



"7 event. The predicted 



photon energy is obtained from the beam energy and 
kinematic information from all particles in the event ex- 
cept the measured photon energy. We observe that the 
excess of data over non-Bhabha MC events is clustered at 
low A£fcy, where the Bhabha events are expected to ap- 
pear. As we progressively loosen the electron veto on the 
tag-side track, the excess in the number of data events 
over the non-Bhabha MC background grows in the region 
with small AE y , providing further confirmation that the 
Bhabha background is well understood. 
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FIG. 2: rriEC distribution of data (dots), the expected back- 
grounds (histograms) and MC signal (curve with arbitrary 
normalization) for \AE - (AE}\ < 2a(AE). 



We cross check the Bhabha contamination in the 
| AS — (AE) | < 2a(AE) region from a data sample with- 
out a tag-side electron veto, by removing the E/p require- 
ment on the tag-side. To estimate the Bhabha contami- 
nation surviving our final event selection, which includes 
a cut of tag-side E/p < 0.5, we use the data in the ad- 
jacent Bhabha-dominated E/p region, 0.5 < E/p < 1.2. 
We extrapolate the rate from the 0.5 < E/p < 1.2 region 
to the E/p < 0.5 region, using a high statistics and high 
purity Bhabha control sample obtained by reversing the 
requirement on (ml/1.8 GeV 2 /c 4 ) - In (2 x pj niss /*/£)/ 2.0 
given above. We estimate the residual Bhabha contam- 
ination in our final selection by multiplying the number 
of events in the 0.5 < E/p < 1.2 region of the no tag-side 
electron veto sample by the ratio of the number of events 
in the Bhabha control sample in the E/p < 0.5 region to 
that in the 0.5 < E/p < 1.2 region. This method gives an 
estimate of 10.3±1.1 Bhabha events inside the ±2a(AE) 
band once the tag-side electron veto is applied. 

In this band, we expect 12.9±2.5 events from the non- 
Bhabha MC backgrounds, thus obtaining a total back- 
ground estimate of 23.2±2.7 events. This compares well 
with the 25 events observed inside the ±2a(AE) band in 
the data. We also find good agreement between the ob- 
served and expected number of events separately for the 
sub-samples with one and three tracks on the tag-side. 

For the final background estimate we use the mEC dis- 
tribution of data events inside the ±2a(AE) band, as 
shown in Figure [21 along with the signal shape included 
for illustrative purposes. The backgrounds from data in- 
side the ±2a(AE) band with \m,EC — m r\ > 3ct(toec) 
are fitted to different orders of polynomials in toec using 
a maximum likelihood approach. A fit with a constant 
probability density function (PDF) yields a total \ 2 of 4.7 
for the 10 bins shown in Figure [21 an d predicts 1.9±0.4 



events inside the final ±2cr(mEc) signal region. Equally 
acceptable goodness of fit is obtained with higher-order 
polynomials. However, the coefficients of the higher or- 
der terms are statistically compatible with zero. The 
background predictions from these PDFs agree with the 
prediction from the constant PDF to within ±0.3 events. 
As these deviations are smaller than the statistical error 
on the prediction from the constant PDF, we conclude 
that the data m-EC distribution is consistent with being 
uniform. 

A cross check using non-Bhabha MC background con- 
tributions combined with residual Bhabha contamination 
estimates obtained from the data is also found to be rea- 
sonably uniform in m^c (Figure[2J and predicts 1.7±0.2 
events inside the ±2cr(mEc) signal box. 

The (5.9±3.7)% difference between data and r-pair 
MC predictions also provides a measure of our ability 
to model the signal- like events in the GSB, since these 
data events have very similar characteristics to the signal, 
both in terms of the trigger response of the experiment as 
well as for the distributions of all the selection variables 
apart from mEC and AE. The systematic error due to 
a particular cut is taken as the product of the marginal 
efficiency of the cut and the relative discrepancy between 
data and MC in the GSB after all other cuts have been 
applied. The contributions from all the different cuts 
added in quadrature yield a 2.3% relative systematic er- 
ror, the only appreciable effect being associated with the 
requirements on m 2 v and p^iss- This approach yields a 
more conservative estimate of the systematic uncertainty 
on the signal efficiency than the more traditional ap- 
proach derived from considering the difference between 
the data and MC prediction for each selection variable, 
which gives a total estimate of 2.0% relative contribution 
from all the cuts. 

The relative systematic uncertainties on the trigger 
efficiency, tracking and photon reconstruction efficien- 
cies, and particle identification are estimated to be 1.4%, 
1.3%, 1.8% and 1.3%, respectively. The requirement that 
the events fall within the ±2cr signal box in rriEC and AE 
contributes a 4.4% systematic error associated with the 
scale and resolution uncertainties of these variables and 
a small contribution from the beam energy uncertainty. 
As we use 1.3 million MC signal events, the contribu- 
tion to the uncertainty arising from signal MC statistics 
is negligible. Adding the contributions of the individual 
terms in quadrature with an additional 2.3% normaliza- 
tion error on the product La TT gives a 6.2% total relative 
systematic uncertainty on ha TT e in the signal box, where 
the efficiency is e — (4.7±0.3)%. We note that our final 
limit on the branching ratio is insensitive to the system- 
atic uncertainty as long as this uncertainty is below 10%. 

We find one event in the signal box for an expected 
background of 1.9±0.4 events. Because of the low back- 
ground levels, we do not fit for a signal in the wec dis- 
tribution as is done in our recent search for r — > /x 7 
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FIG. 3: 90% C.L. upper limits on Mf^/Mf^ for Bfr* 
e ± 7) < 1.1 x 10~ 7 with tan/3 = 10, 20 and 40. 



Rather, we set an upper limit employing the same 
technique used in our sear ch for t± -> [3 where 

the background levels were also small. A 90% C.L. up- 
per limit on the branching ratio is calculated according 
to Bfj° L = Nl r L /{2eha TT ) 1 where N$° L is the 90% C.L. up- 
per limit with one event observed when 1.9±0.4 events 
are expected. The limit is calculated including all uncer- 
tainties using the technique of Cousins and Highland |io| 
following the implementation of Barlow 21]. At 90% 
C.L. this procedure gives an upper limit of S(r ± — > e ± 7) 
<l.lxlO -7 [2^ ] . This represents a more than three-fold 
reduction in the upper limit as reported in [jj. 

As an example of how this result constrains theories 
beyond the SM, we set bounds on the ratio of the first 
and the third generation element to the first generation 
diagonal element (-^f 13 /-^f n ) °f the left-handed slep- 
ton mass matrix based on predictions from a minimal 
supergravity model [U [24| . Figure |3| shows the upper 
limits on M~ /M~ as a function of the ratio of the 
Higgs vacuum expectation values (tan p) and the univer- 
sal scalar mass (mo), which, for simplicity, is set equal to 
the universal gaugino mass. 
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A search for the non-conservation of lepton flavor in the decay r — * e ± j has been performed with 2.07 x 10 s 
e + e~ — * t + t - events collected by the BABAR detector at the PEP-II storage ring at a center-of-mass energy near 
10.58 GeV. We find no evidence for a signal and set an upper limit on the branching ratio of B{t ± — > e ± 7) < 1.1 x 10~ 7 
at 90% confidence level. 



